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How can we quantify the entanglement in a quantum state, if only the expectation value of 
a single observable is given? This question is of great interest for the analysis of entanglement in 
experiments, since in many multiparticle experiments the state is not completely known. We present 
several results concerning this problem by considering the estimation of entanglement measures via 
Legendre transforms. First, we present a simple algorithm for the estimation of the concurrence and 
extensions thereof. Second, we derive an analytical approach to estimate the geometric measure of 
entanglement, if the diagonal elements of the quantum state in a certain basis are known. Finally, 
we compare our bounds with exact values and other estimation methods for entanglement measures. 

PACS numbers: 03.65.-w, 03.65.Ud, 03.67.-a 



I. INTRODUCTION 

Entanglement is a key phenomenon in quantum infor- 
mation science and the quantification of entanglement is 
one of the major problems in the field. For this quantifi- 
cation many entanglement measures have been proposed 
P, [3, i, i, i, i, m, S, il. However, a central problem 
in most of these proposals is the actual calculation of a 
given measure: entanglement measures are typically de- 
fined via optimization procedures, which may consist of a 
maximization over certain protocols or the minimization 
over all decompositions of a state into pure states. For 
some remarkable cases it happens that such minimiza- 
tions can be performed analytically [13, [HI, [H, [3 1 how- 
ever, in the general case these problems are not solved. 
Therefore, to take a realistic point of view, one can try 
to estimate entanglement measures, and many proposals 
for the estimation of ent ang lement measures have been 
presented [H, [H, [H, [13, [ill • 

In an experimental setting, the situation becomes 
even more complicated: since quantum state tomog- 
raphy for multiparticle systems requires an exponen- 
tially increasing effort, the state is often not completely 
known. Typically, one measures so-called entanglement 
witnesses, special observables, for which a negative ex- 
pectation value signals the presence of entanglement 
In, [13, [m, [13, [11,111, [H]. In this sense, entanglement 
witnesses allow to detect entanglement, but the ques- 
tion arises, whether they also allow to quantify entan- 
glement. This question has been addressed from several 
perspectives [11 [13, [H, [13, [lOl and in Refs. [sil, [H a 
general recipe for this problem was found. There it has 
been shown how one can derive the optimal lower bound 
on a generic entanglement measure from the expectation 
value of a witness or another observable. The estimate 
uses Legendre transforms to give lower bounds on the 
convex entanglement measure, and the main task in this 
scheme is to compute the Legendre transform of a given 



entanglement measure. 

In this paper, we extend this method into several di- 
rections. First, we derive a simple algorithm for the cal- 
culation of the Legendre transform for the concurrence 
d] and extensions thereof. Then, we present analytical 
results for the Legendre transform for certain witnesses 
for the geometric measure of entanglement [3. Finally, 
we discuss examples and compare our results to other 
methods for entanglement estimation. But before pre- 
senting the new results, let us shortly review the method 
presented in Refs. [3ll.[3^. 

II. THE METHOD 

Let us consider the following situation: in an exper- 
iment, an entanglement witness W has been measured 
and the mean value (W) = Tr{Q'W) = w has been found. 
The task is now to derive from this single expectation 
value a quantitative statement about the entanglement 
present in the quantum state. In our case, we aim at 
providing a lower bound on the entanglement inherent in 
the state g. That is, we are looking for statements like 

Tr{Q-W)=w => E{g)>f{w), (1) 

where E{g) denotes an arbitrary convex and continuous 
entanglement measure. We do not specify it at this point 
further. Naturally, we aim to derive an optimal bound 
/(w) and an estimate is optimal, if there is a state go 
with Tr{goW) = w and E{gQ) = f{w). 

In order to derive such lower bounds, let us consider 
the so-called Legendre transform of E for the witness W, 
defined via the maximization 

E{W) = sup{Tr(>Vg) - E{g)}. (2) 

Q 

As this is defined as the maximum over all g, we have for 
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FIG. 1: Geometrical interpretation of the Legendre transform. 
The dotted line is a general affine lower bound, and the dashed 
line is an optimized bound with the minimal /i. By varying 
the slope, one finds an affine lower bound, which is tight for 
a given wi (solid line). See text for details. 

any fixed g that E{W) > Tr{Wg) - E{g), hence 

E{g) > Tr{Wg) - E{W), (3) 

which is known as Fenchcl's inequality or Young's in- 
equaUty. The point is that the first term on the right 
hand side are the given measurement data, while the 
second term can be computed. Therefore, a measurable 
bound on E{g) has been obtained. 

In order to improve this bound, note that knowing the 
data Tr{gW) = w is, of course, equivalent to knowing 
Tr{gXW) ~ Xw for any A. Therefore, we can optimize 
over all A and obtain 

E{g) > sup{XTr{Wg) - E{XW)}. (4) 

A 

This is a better bound than Eq. ^ and, as we will see, 
already the optimal bound in Eq. 

For our later discussion, it is important to note that 
this estimation has a clear geometrical meaning [see 
Fig. (1)]. As E{g) is convex, the minimal E{g) com- 
patible with Tr(gW) = w, denoted by f{w) is convex, 
too HI]. Let us consider a generic afSne lower bound 
g{w) = Xw — /i on f{w), i.e., /(w) > Xw — fi. We have 
II > Xw — f{w), and in order to make the bound for a 
fixed A as good as possible, we have to choose /z as small 
as possible. This leads to fiopt — sup^{Xw— f{w)} which 
is exactly the optimization in Eq. 

For a fixed slope A we obtain by this method an afhne 
bound (characterized by ^lopt) which is already optimal 
for a certain value wq, as it touches f{w) in this point. 
For any other mean value w it delivers a valid, but not 
necessarily optimal bound. To obtain the optimal bound 
for any given wi we have to vary the slope A. This corre- 
sponds to the optimization in Eq. Q. Since / is convex, 
we obtain for each w the tight linear bound, showing that 



this optimization procedure gives indeed the best possi- 
ble bounds in Eq. ([1]). 

Three remarks are in order at this point. First, it was 
not needed that W is an entanglement witness, we may 
consider an arbitrary observable instead. Second, we can 
also consider a set of observables VV = {Wi, yy„} 
at the same time. We just have to introduce a vector 
A = {Ai,...,A„} and replace AW by J^k^kyVk, then all 
formulas remain valid. Finally, also for a non-convex 
E{p) the method delivers valid bounds, however, then it 
is not guaranteed that the bounds are the optimal ones. 

In any case, the method relies vitally on the ability to 
compute the Legendre transform in Eq. The diffi- 
culty of this task clearly depends on the witness W and 
on the measure E{g) chosen. At first sight, the task may 
seem hopeless, as the calculation of E{g) for mixed states 
is for many measures already impossible. 

For instance, a large class of entanglement measures is 
defined via the convex roof construction. For that, one 
first defines the measure E{\tJj)) for pure states, and then 
defines for mixed states 

inf Vpfcii;(|0fc)), (5) 

k 

where the infimum is taken over all possible decom- 
positions of g, i.e., over all pk and \cj)k) with g = 
J2kP''\4'k){4'k\- Clearly, this infimum is very difficult to 
compute. 

However, for the computation of the Legendre trans- 
form, this is not relevant: as one can easily prove (see 
Ref. [3l| for details), for convex roof measures the maxi- 
mization has to run only over pure states 

^(W)=sup{(V'|W|V)-i?(|V'))}, (6) 
W 

which simplifies the calculation significantly. In fact, this 
shows that for convex roof measures, which are by con- 
struction rather difficult to compute, the Legendre trans- 
form is rather simple to compute. 

In Ref. [m we have considered the entanglement of 
formation and the geometric measure of entanglement 
as entanglement measures, which are both convex roof 
measures. We have provided simple algorithms for the 
calculation of the Legendre transform, and for special 
witnesses, we calculated the Legendre transform of the 
geometric measure also analytically. In this paper, we 
will consider the concurrence and extensions thereof, and 
we will also present new analytical results for the geomet- 
ric measure. 

Finally, it should be noted that the optimization in 
Eq. ^ can be completely skipped. Any A delivers al- 
ready a valid bound. However, the optimal A can easily 
be found numerically. 
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III. THE CONCURRENCE 



So we make the ansatz 



As a first entanglement measure, let us discuss the con- 
currence. This quantity is defined for pure states as [3| 



Ecm = ^2[l -Tr{g\)], 



(7) 



where qa the reduced state of {ip) for Alice. For mixed 
states this definition is extended via the convex roof con- 
struction in Eq. ([5]). For the special case of two qubits, 
the concurrence is a monotonous function of the entan- 
glement of formation, moreover, the minimization in the 
convex roof construction can be explicitly performed . 
In the general case, it is not so directly connected to the 
entanglement of formation. The concurrence is an entan- 
glement monotone [IHl, but it does not fulfill all desired 
axioms for an entanglement measure, e.g. it is not addi- 
tive. 

Here, we want to calculate the Legendre transform of 
Ec for a generic witness W. Similar as it has been done in 
Ref. [sH for the entanglement of formation, we will con- 
struct an iterative algorithm for the optimization, which 
converges to the maximum as a fixed point. 

First, as the concurrence is defined via the convex roof, 
it suffices in Eq. ^ to optimize over pure states only. 
Then, using the fact that for x E [0; 1] 



Qe[0;i] 2a 2 
we can rewrite the Legendre transform as 



Ec{y^) = sup sup{(V'|>V|?/') 



V2a 



-— }• 



(8) 



(9) 



The idea is to write this maximization as an iteration that 
optimizes a and {ip) in turn. For a fixed we perform 
the maximization over a analytically and similarly we can 
find the optimal [ip) for a fixed a. Concerning the first 
step, note that for a fixed 1-0) the optimal a is simply 
given by 



a 



(10) 



Concerning the second step, if a is fixed, we have essen- 
tially to solve an optimization problem like 

sup{(V|W|V')-[l-^r(^^^)]}, (11) 

w 

where W is proportional to the original witness W. 

Here, the second term is nothing but the q-entropy, 
investigated in detail by Havrda, Charvat, Daroczy and 
Tsallis fllllllli,!!!. 



Sqig) = 



1 - Tr{Q'i) 



(12) 



for the case q = 2. Therefore, we can try to write Sq{gA) 
as an infimum via the Gibbs principle, similarly as it has 
been done for the von Neumann entropy in Ref. [sij . 



S2{g) = M{Tr{QH) - F^iH)} 
F2{H) = mi{Tr{QH) - S2{q)} 

Q 



(13) 
(14) 



For the case of the von Neumann entropy these formulas 
just express the Gibbs variational principle, where F is 
the free energy, and the inverse temperature was set to 

We are mainly interested in the first minimization and 
have to compute F2 and the H, where the minimum is 
attained. The point is that the second minimization has 
been solved already by Gucrbcroff and Raggio [i^ and 
the unique thermal state minimizing Eq. (jl4p for an arbi- 
trary Hamiltonian H and consequently F2 is known. For 
the first minimization, it remains to find the Hamiltonian 
for which the given state is the thermal state. 

In order to do this in practice, let us first recall the 
results of Ref. [i^l- For the case q = 2 and /3 > the 
results of this reference state the following: 

Let H he a Hamiltonian with ground state energy S- , 
ground state degeneracy g- and let the energy of the first 
excited state be el . Define 



13' 



> and r 



1 



> 0, 



and the monotonically increasing function p(t) : [0, 
[0,(3+] as 



(15) 



m = 



2t 



Tr([l -i(i?-£_l)]e 



(16) 



where [X]^ denotes the positive part of the operator X. 
Let T be the inverse function to P{t). 

Then, for /? < (3+ the unique thermal state is given by 



Q = N{[l^T{P){H-e^l)], 



(17) 



where M denotes the normalization, and for (3 > (3+ the 
thermal state is just the normalized projector onto the 
eigenspace corresponding to the lowest energy e_. From 
this, it is clear that g and H arc diagonal in the same 
basis. 

Coming back to the minimizations in Eqs. (fT3l [T4| let 
us assume that a density matrix g with decreasing eigen- 
values \i (i = 1, ...,iV) is given, and the task is to com- 
pute the corresponding Hamiltonian H with increasing 
eigenvalues Ei . Without loosing generality, we can choose 
El = 0. 

Let us first consider the case that g is non degen- 
erate and has full rank. Then, Eq. p7|) implies that 
the eigenvalues Ei of H have to fulfill 1 — T{(3)Ei > 
for all i. Since we are considering the case /3 = 1, 
we have for tq = t(/3 = 1) due to Eq. that 
To = [Tr(l - ToH)]/2. From Eq. ^ is follows first 
that Ai = TV = 1/Tr{l - toH) = 1/(2to), and then 
A, = Ai(l - TnE,) = Ai[l - E,/{2Xi)], hence 



E, = 2(Ai - A,). 



(18) 
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For this Hamiltonian, wc have also 

F2(if) = 2Ai-^A2-l. (19) 

i 

As a remark, first note that this solution fulfills 1 — 
T{f3)Ei = 1 — (Ai — Ai)/Ai > 0, as requested at the bc- 

I 



If and H are fixed, we can compute the optimal a 
as shown in the beginning. If \ip) and a are fixed, wc 
can determine the optimal H as above. Finally, if H and 
a are fixed we choose \ip) as the eigenvector correspond- 
ing to the maximal eigenvalue of W — {H (g) l)/{\/2a). 
Therefore, we have an iterative optimization, which de- 
livers an monotonously increasing sequence of values for 
the Legendre transform, with the actual value as a fixed 
point. 

This algorithm can be implemented with few lines of 
code, examples will be discussed in Section V. It should 
be noted that in principle we cannot prove that the al- 
gorithm converges always to the global optimum, as the 
final fixed point may depend on the starting values of 
\ip), H and a, leading to an overestimation of the entan- 
glement. In practice, however, the convergence behavior 
is quite good naturcd and the algorithm delivers a well 
suited tool for obtaining sharp bounds on the concur- 
rence. 

Finally, let us add that with the same method also 
extensions of the concurrence may be treated. For in- 
stance, one may consider quantities for pure states as 
^qiW) = SqiQA)- With the convex roof extension, these 
are also entaiiglement monotones [ssj . and using the re- 
sults of Ref. [40| , the Legendre transform can be calcu- 
lated in a similar manner as above. 

In addition, there are several multipartite entangle- 
ment monotones, which can be seen as multipartite ex- 
tensions of the concurrence 0, H, [4l|. For instance, for 
A^-qubit states one can define the Meyer Wallach measure 



EMwm^2[l--Y.^rigl)], (21) 
fc=i 



where the gk are all reduced one-qubit states. Clearly, 
due to the similar structure as in Eq. ([7]) the methods 
from above can also be used to compute the Legendre 
transformation for these measures. 



ginning. Second, it delivers /?+ = l/(Ai — A2) > P = 1, 
justifying the ansatz in Eq. ([T7|) . Finally, it is easy to see 
that if g is not of full rank or degenerated, the recipe for 
above works also and delivers a correct solution. 

In summary, we can write the problem of the compu- 
tation of the Legendre transform for the concurrence as: 



(20) 



IV. THE GEOMETRIC MEASURE 

As the second entanglement measure, let us consider 
the geometric measure of entanglement Eg 0. This is 
an entanglement monotone multipartite systems, quanti- 
fying the distance to the separable states. The geometric 
measure is defined for pure states as 

Eaim^l- sup (22) 

|0> = |a)|&)|c)... 

i.e., as one minus the maximal overlap with pure fully 
separable states, and for mixed states via the convex roof 
construction. 

The geometric measure is a multipartite entanglement 
measure, as it is not only a summation over bipartite en- 
tanglement properties. Despite of its abstract definition, 
it has turned our that Eq can be used to quantify the 
distinguishability of multipartite states by local means 
[4^ . As the geometric measure is one of the few measures 
for multipartite systems, which have a reasonable opera- 
tional meaning and are at the same time proved to fulfill 
all of the conditions for entanglement monotones, it has 
been investigated from several perspectives, for instance 
it has been used to study multipartite entanglement in 
condensed matter systems Q- 

In Ref. [3l| the problem of calculating the Legendre 
transform for the geometric measure was already consid- 
ered and the following results were obtained: 

First, an iterative algorithm (in the same spirit as the 
algorithm for the concurrence in the previous section) 
has been derived for calculating the Legendre transform 
of arbitrary witnesses. Second, for the important case 
that the witness is of the form W = al — |x)(xl (or, 
equivalently W = |x)(x|) an analytic formula of the Leg- 
endre transform has been derived, reading, 

Ea{rW) = 

f ra for r > 0, 

\ [y/{l - r)'^ + irEGilx)) + 2ar - r - 1] /2 for r < 0. 

(23) 

Here, we want to generalize this result by determining 
analytically Eg{W) for the case that W is diagonal in 



i?c(W) = sup sup supf(V|W|V)--^[(V'|(i?®l)|^)-i^2(H)] -^|. 
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some special basis, e.g. the GHZ-state basis. We first 
consider the special case of the GHZ state basis, then 
the formula for the general case can directly be written 
down. 

So let us assume that 



W = Y,^^\GHZ,){GHZ, 



(24) 



i=l 



is the witness for an N qubit system, where \GHZi) = 
± |y(i)...y(^)))/V2 with € {0,1} 

and a;*^*^ ^ j/*^*^ is the GHZ-state basis. Without loos- 
ing generality we assume that the Ai are decreasingly 
ordered, i.e., Ai > A2 > ... > A2iv, but not necessarily 
positive. 

The Legendre transform is given by 



sup sup 

\4>) = \a)\b)\c). 



\[w- 



sup 

) = \a)\b)\c). 



1, (25) 



where \\X\\ denotes the maximal eigenvalue of the oper- 
ator X. In order to compute this maximal eigenvalue, we 
write the operator [W-l- 10) in the GHZ basis. W is di- 
agonal there and since the maximal overlap between the 
fully separable state \4>) and any of the GHZ states is 1/2 
(i.e., the geometric measure for GHZ states is 1/2 [3|), 
the matrix representation of |0)((/>| has matrix elements 
with absolute values not larger than 1/2. 

Our claim is now that the optimal choice of \4>) is to 
take 10) (01 as a 2 X 2 matrix with all entries 1/2 and 
acting on the two-dimensional space corresponding to the 
largest eigenvalues Ai and A2. To prove that this choice is 
really optimal, we show that the above mentioned choice 
is also optimal if we consider the more general class of all 
10), which have an overlap smaller or equal 1/2 with all 
the GHZ states, but which are not necessarily product 
states. 

We prove it by contradiction. Let us assume a different 
optimal solution |0) and a corresponding eigenvector to 
the maximal eigenvalue IV')- The vectors can be written 
as 10) = EfeafclGiJ^fc) and \^) = Y.kh\GHZk). We 
assume without loosing gener ality that < jcup < 1/2 
and < I asp < 1/2. The fimction to maximize is given 
by 



Earn = \Y.alh\^ + |/3fcpAfe - 1 



(26) 



Since we are interested in the maximum, we can, without 
restricting generality assume that au and are real and 
positive. The interesting terms for the discussion are 



(27) 



X is a scalar product, which is maximal if the vec- 
tors (ai,Q;3) and {Pi, (is) are parallel. For given val- 
ues of (/3i,/?3) this may be prohibited by the constraint 



< 1/2. Then, however, it is clearly optimal to take an 
(ai^a^) at the border of the domain, which has for one 
ai = 1/2, leading to a contradiction to the assumption on 
the form of |0). Otherwise, we can choose the two vectors 
parallel, and (01,03) is not at the border. Then, how- 
ever, we can enlarge Y by increasing (3i in (/?i,/33) (and, 
simultaneously ai in (ai, a^) in order to keep the vectors 
parallel). This leads to another contradiction concerning 
the optimality of \tp). 

Note that the class of the considered |0) is strictly 
larger than the class of product vectors, since not for any 
pair of GHZ states \GHZi) and \GHZ2} we can find a 
product vector |0), such that |0) has an overlap of 1/2 
with both of the \GHZi). However, we obtain an upper 
bound on the Legendre transform from this ansatz, which 
can be used for a valid lower bound on the entanglement 
measure. Further, one can check whether this bound is 
tight by direct inspection of \GHZi) and \GHZ2) after- 
wards. 

Having shown that the simplest choice of |0) is optimal, 
the calculation of the Legendre transform reduces to a 
calculation of eigenvalues of a 2 x 2-matrix, and we have: 



Ai + 1/2 
1/2 

Ai + A2 - 1 



1 



1/2 
A2 + 1/2 

+ ^V(Ai-A2)2-L (28) 

The generalization to other states besides GHZ states 
is straightforward: if the overlap is bounded by some 
other number (e.g. 1/4 for four-qubit cluster states), we 
only have to calculate the eigenvalues of some larger ma- 
trix (e.g. a 4 X 4-matrix for four-qubit cluster states), 
in order to derive an analytical upper bound on £'g(W). 
We can summarize: 

Observation. Let W = X]i=i •^ilV'j)('0i| be an op- 
erator, where for all eigenvectors 1-0^) the overlap with 
fully separable states is bounded by l/k,k G N. Then 
the Legendre transform is bounded by 



EciW) <\\[X] 



(29) 



where X is a k x fc-matrix with the entries Ai -I- 1/fc on 
the diagonal, and offdiagonal entries l/k. The question 
whether this bound is the exact value, can be decided by 
direct inspection of the 

This Observation allows for a simple calculation of a 
lower bound on Eq if the fidelities of the basis states 
are known. First, the estimation is much simpler 
and faster compared with the iteration algorithm for ar- 
bitrary witnesses, since the optimization runs only over 
the k largest eigenvalues A^ of the possible witnesses 
W ~ J2i •^i|'0i)(V'i|- For the iteration algorithm, it would 
be necessary to consider all witnesses W = J2i ^ili^i) ii^i I j 
which amounts to a variation over 2^ parameters A^. Sec- 
ond, the bounds on Eg may become significantly better, 
compared with the estimation from a single fidelity, ac- 
cording to Eq. 

For example, the |V'i) may be a graph state basis, where 
the fidelities have been determined from the expectation 
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values of the stabilizer operators [23|. In Section V we 
will discuss an example for four-qubit states. 



V. EXAMPLES 

In this section, we present several examples for the 
presented method and compare it with other estimation 
methods as well as with exact values of the entanglement 
measures. A Mathematica file with the used algorithms 
for the calculation of the Legendre transforms is available 
from the authors. 



A. Concurrence for isotropic states 

As a first example, let us consider isotropic states in a 
N X N system, defined by 

g{F) = ^j^{l-\cj^){cb\)+F\cj^){cbl (30) 

which are a convex combination of a maximally en- 
tangled state 10) = ^^ilii)/ VN and the totally mixed 
state. The parameter F encodes the fidelity of i.e., 
F ~ ((j)\p(F)\cj)). For these states, the concurrence is 
known to be [l3 



Cig) 



N-V n' 



(31) 



In order to test our methods, we consider the standard 
witness for states of the for states of the from in Eq. (pO|) . 
namely 



1 

N 



I0>('/'I, 



(32) 



and estimate from its expectation value the concurrence, 
using our algorithm. The results for the case = 3 are 
shown in Fig. (2). It turns out that for this case, our 
lower bounds are sharp and reproduce the exact value of 
the concurrence. 



B. Comparison with other estimation methods 

Let us compare the presented estimation method with 
other methods of estimating entanglement measures. For 
this aim, we consider two-qubit states of the form 



g{p) =p|V')(V'| + (1 -p)cr, 



(33) 



where {ip) is a pure entangled state and a is unknown and 
random separable noise. As entanglement measure we 
consider the entanglement of formation. This is, for the 
case of two qubits, equivalent to the concurrence, since it 
is a monotonous function of it. We consider six different 
methods for estimating the entanglement of formation: 



1.2 



0.8 



0.6 



0.4 



0.2 



Exact value — 
Lower bound □ 




















































- 































0.3 



0.4 



0.5 



0.6 0.7 
F 



0.9 



FIG. 2: Exact value of the concurrence and a lower bound 
from the witness in Eq. (|32|l and the algorithm in Sec. IV for 
isotropic states. See text for details. 



1. For the two-qubit case we can exactly calculate the 
entanglement of formation due to the formula of 
Wootters ■ Clearly, since this requires complete 
knowledge of the density matrix, for an experimen- 
tal implementation state tomography is needed. 

2. In Refs. [13, a method to estimate the entan- 
glement of formation or the concurrence from the 
separability criterion of the positivity of the partial 
transpose (PPT) or the computable cross norm or 
realignment criterion (CCNR) was presented. Ex- 
perimentally, this approach requires again state to- 
mography. 

3. We can also take the witness which is proper for the 
states of the form e(p) = p\ip){ip\ + {l-p)l/4. This 
witness might not be the optimal one for the state 
under investigation, since the noise is not known 
and in general white. However, we can use the 
Legendre transform with the algorithm of Ref. [3l| 
to estimate the entanglement of formation from it. 
Equivalently, we can use the algorithm of Section 
III to estimate the concurrence. Experimentally, 
this method does not require state tomography, 
only three local measurements are needed for the 
measurement of the witness (4^ . 

4. The witness in the third method is of the form 

W-|0)(<^r-. (34) 

Clearly, the mean value (W) = rr(goxpVV) = 
{(j)\g'^^p\(j)) can be used to derive a lower bound on 
the negativity of the partial transpose g^^p- Then, 
the second method [3, may be used to estimate 
the entanglement of formation. This approach does 
not require state tomography. 

5. In Ref. lower bounds on the concurrence from 
measurements on two copies of the state g have 
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FIG. 3: Comparison of different methods to estimate the en- 
tanglement of formation for noisy two-qubit states. EXACT 
denotes the exact calculation according to Ref. [l^l; CAFl the 
second method, estimating the entanglement form the PPT 
criterion [ij, Oi WIT the estimation from the witness as 
presented in this paper and in Ref. ^31]; CAF2 the fourth 
method as a combination of the witness with CAFl: 2COP 
the fifth method, using measurements on two copies and 
RWIT the sixth method [H. See text for details. 



Method 


CAFl 


WIT 


CAF2 


2C0P 


RWIT 


# of measurements 


9 


3 


3 


9 


3 


?7 for p G [0.8; 1.0] 


96.0 % 


94.8 % 


94.8 % 


70.1 % 


63.0 % 


7? for p G [0.6; 0.8] 


86.0 % 


72.0 % 


72.0 % 


5.9 % 


7.0 % 



TABLE I: Comparison of the different estimation methods. 
For each method, the required number of local measurement 
settings on a single copy and the efficiency rj is given. Here, 
the efficiency is defined as the ratio between the estimated 
value of the entanglement of formation and the actual one, 
for two different regions of the noise parameter p. 



which are very close to the exact result. The second 
method is the best bound, however, it requires complete 
knowledge of the state. The third method is by con- 
struction better than the fourth method (as the Legen- 
dre transform delivers by construction the best possible 
bounds from a given witness). In this example, how- 
ever, they are practically equivalent. The fifth and sixth 
method deliver good results if g is close to a pure state. 



C. Geometric measure for four-qubit states 



been derived. For the case of two qubits we can 
use them to bound also the entanglement of for- 
mation. A measurement of a single observable on 
two copies can always be expressed as a function of 
mean values of local obscrvables on a single copy. 
In this case, however, for such an implementation 
effectively state tomography is needed. 

6. In Ref. [sOj a method for estimating the concur- 
rence from special entanglement witnesses has been 
derived. Namely, it was shown that 

i^c(.)>-^(>V,), (35) 

where W4, = 2[1^ (g) Tr^d^) (0|) - |0)(0|] is a wit- 
ness belonging to the reduction criterion. We use 
this method with the witness for the state 
With this choice, the method automatically repro- 
duces the exact value for the case p = 1 and it can 
be seen that Eq. ([35]) is nothing but the Legendre 
transform for a special choice of the slope A. Ex- 
perimentally, this method would also require three 
local measurements [13, El] ■ 

As an example, we considered in Eq. ([33|l the pure state 
— (4 1 00) + |ll))/\/T7 and the state a was randomly 
(in the Hilbert-Schmidt measure) chosen from the set of 
separable states [1^. For a fixed a we calculated all the 
above mentioned values depending on the noise level p 
and finally averaged over hundred realizations of a. The 
results are shown in Fig. (3) and Table 1. 

One can clearly see that the methods two, three and 
four result in bounds on the entanglement of formation, 



As a third example we discuss the geometric measure 
of entanglement. In order to demonstrate the method in 
Section IV, we consider an experimental situation similar 
to the one in Ref. [2^. In this experiment, a four-photon 
cluster state 

\CL) = i(|0000) + jOOll) + jllOO) - |1111)) (36) 

has been prepared using parametric down conversion and 
a controlled phase gate. The fidelity of the target state 
was then determined by the measurement of all stabilizer 
operators. These operators are given by the local obscrv- 
ables 5"! = a^a^U, 5*2 = axCTxCrz'^, = Ict^ct^CTj;, and 
54 = IIctzCt^, and products of these observables. The 
cluster state is an eigenstate (with eigenvalue -1-1) of all 
these 2^^ = 16 obscrvables, and from their expectation 
values the fidelity can be determined [U, [2^ . 

Using the fact that the maximal overlap of the cluster 
state with fully separable states equals 1/4 (i.e., the geo- 
metric measure is 3/4 [i^l) Eq. (|^ can be used to bound 
Eq from this fidelity. There are, however, also other 
common eigenstates of the Si with different eigenvalues. 
These states are orthogonal to the cluster state and form 
the so-called cluster state basis. All states in this basis 
share the same entanglement properties and their fideli- 
ties can also be determined from the mean values of the 
Si. In fact, the knowledge of all {Si) is equivalent to the 
knowledge of all fidelities. 

In order to investigate how the information about the 
fidelities of all states in the cluster state basis can be 
used for the estimation of entanglement, we consider the 
simple case that only three fidelities are larger than zero, 
Fi . F2 and F3 = 1 — Fi — F2 . For a given triple of fidelities 
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FIG. 4: Comparison between two estimation methods for 
the geometric measure of entanglement. The lower curve is 
the analytical bound from the maximal fidelity according to 
Ref. [sJl or Eq. ((23}, the upper curve is the analytical bound 
from Section IV, taking all fidelities into account. See text 
for details. 

we may first consider the maximal fidelity and then use 
Eq. (|23l) to obtain a lower bound on Eq. Alternatively, 
we can use the methods of Section IV and consider all 
three fidelities at the same time. In practice, this gives a 
lower bound on Eq by the optimization problem 

3 

Eg> sup {^AfcFfe-||[X]||+l}, (37) 

Ai,...,A4 j.^^ 



where X is defined as in Eq. ([29]) for fc = 4. Any set of Xi 
delivers already a valid lower bound, and the optimum 
over all Xi is easily found. 

The results are plotted in Fig. (4). One can clearly see, 
that taking all fidelities into account, improves the lower 
bounds significantly. 



VI. CONCLUSION 



In conclusion, we have investigated how entanglement 
measures can be estimated from incomplete experimen- 
tal data. We have shown that the method of Legendre 
transforms can successfully be applied to the concurrence 
and extensions thereof. Furthermore, we have presented 
an analytical way to estimate the geometric measure if 
the fidelities of certain basis states are known. Extending 
the presented methods to other entanglement measures 
is an interesting task for further study. 

We thank J. Eisert, N. Kiesel, A. Osterloh, J. Siewert 
and W. Wieczorek for discussions. This work has been 
supported by the FWF, the DFG, and the EU (OLAQUI, 
PROSECCO, QUPRODIS, QICS, SCALA). 
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